Kaon Condensation in Dense Matter 



J. Carlson 1 , H. Heiselberg 2 and V. R. Pandharipande 3 

1) Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
2) NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0., Denmark 
3) Dept. of Physics, Univ. of Illinois at Urbana- Champaign, 1110 West Green St., Urbana, Illinois 61801-3080, USA 

(February 9, 2008) 



o 

Q 



> 

m 
o 

ON 



o 

=5 



X 



The kaon energy in neutron matter is calculated analyti- 
cally with the Klein-Gordon equation, by making a Wigner- 
Seitz cell approximation and employing a K~ N square well 
potential. The transition from the low density Lenz potential, 
proportional to scattering length, to the high density Hartree 
potential is found to begin at fairly low densities. Exact non- 
relativistic calculations of the kaon energy in a simple cubic 
crystal of neutrons are used to test the Wigner-Seitz and the 
Ericson-Ericson approximation methods. All the calculations 
indicate that by ~ 4 times nuclear matter density the Hartree 
limit is reached, and as the Hartree potential is less attrac- 
tive, density for kaon condensation appears to higher than 
previously estimated. Effects of a hypothetical repulsive core 
in the K~ N potential are also studied. 



Kaon condensation in dense matter was suggested by 
Kaplan and Nelson M, and has been discussed in many 
recent publications Due to the attraction be- 

tween K~ and nucleons its energy decreases with in- 
creasing density, and eventually if it drops below the 
electron chemical potential in neutron star matter in (3- 
equilibrium, a Bose condensate of K ~ will appear. 

The kaon-nucleon interaction in vacuo have been de- 
scribed by Brown, Lee, Rho, and Thorsson Q by an ef- 
fective Lagrangian based on chiral perturbation theory 
The K~n interaction is well described and is fortunately 
not affected much by resonances as is the K~p interac- 
tion. Kaiser, Rho, Waas and Weise [EJ have used en- 
ergy dependent KN amplitudes calculated in a coupled- 
channel scheme starting from the chiral SU(3) effective 
Lagrangian. They correct for correlation effects in a nu- 
clear medium and find that K~'b condense at densities 
above ~ 4po, where po — 0-16 fm -3 is normal nuclear 
matter density. This is to be compared to the central 
density of ~ 4po for a neutron star of mass 1.4M Q accord- 
ing to the estimates of Wiringa, Fiks and Fabrocini || 
using realistic models of nuclear forces. The condensate 
could change the structure and affect maximum masses 
and cooling rates of massive neutron stars. 

In this letter we calculate the kaon energy in neu- 
tron matter using the Wigner-Seitz approximation for 
the Klein- Gordon equation of kaons in neutron matter. 
Our formulation is exact in both, the low density and the 
high density limits. 

We assume that the kaon-nucleon interaction is via the 
Weinberg- Tomozawa vector potential V(r). In the anal- 



ysis of Ref. H the K + N interaction was also found to be 
dominated by u) and p vector mesons. The energy of the 
kaon-nucleon center-of-mass system with respect to the 
nucleon mass is then 
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where mjv = 939.5 MeV is the neutron mass, mx — 
494 MeV the kaon mass, and k is the kaon momentum 
(we use units in which h and c are unity) in center-of-mass 
frame. We have included the recoil kinetic energy of the 
nucleon assuming that terms of order /c 4 /8m 3 v and higher 
can be neglected. For a relativistic description of the 
kaon in a vector potential we employ the following recoil 
corrected Klein-Gordon (RCKG) equation obtained by 
quantizing Eq. (Ill) (k = — iV) 



(co-V(r)f + 



raw + oj — V(r) , 



m N 



= 0. (2) 



The shape of the kaon-nucleon potential is not known. In 
most of our work we approximate it with a square well: 



V(r) = ~V Q{R - r) 



(3) 



Yukawa and repulsive core shaped potentials are less fa- 
vorable for kaon condensation. The range of the inter- 
action R and the potential depth Vq are related through 
the s-wave scattering length: 



a = R- 



tan(Koi?) 



(4) 



where Kq — (2itikVq + V 2 )iriN /{tun + rax + Vq). If 
V « m K, this reduces to the nonrelativistic result Kq — 
2mflVo, where m# = rnxmN /(mx + tun) is the kaon- 
nucleon reduced mass. The K~n scattering length is 
negative corresponding to a positive Vo- 

The kaon-nucleon scattering lengths are, to leading or- 
der in chiral meson-baryon perturbation theory, given by 
the Weinberg- Tomozawa vector term 



a>K±n = a K ± p /2 = ± 



m R 
Svr/2 



±0.31fm, 



(5) 



where / ~ 90 MeV is the pion decay constant. Scalar 
and other higher order terms have been estimated and 
we shall use the value a K - n = —0.41 fm from Qj. This 
estimate does not include the T,tt decay channels which 
account for the complex parts of the K~N scattering 
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lengths |?]|| . The empirical scattering lengths extracted 
from scattering measurements as well as kaonic atoms are 
a K - n = (—0.37 — i0.57) fm and a K - p — (0.67 — i0.63) 
fm. 

Effective Lagrangians predict a^- p ~ —0.82 fm, in 
perturbation theory. However, it becomes positive due to 
the presence of the nonperturbative A(1405) resonance. 
Interestingly, we can estimate the range of our square 
well potential assuming that the A(1405) is a K~p bound 
state. As in || we assume the K ~p potential to be twice 
as strong as the K~n potential given by Eq. (Q). In 
order to form a K~p bound state with binding energy of 
m p +m,K — mA(1405) = 27 MeV, and havea^-„ = —0.41 
fm we need R ~ 0.7 fm. For this reason we present results 
for R = 0.7 fm; however our main conclusions are valid 
for reasonable values of R. Using a — —0.41 fm and 
R = 0.7 fm we obtain V =122 MeV with RCKG and 
=126 MeV with the non-relativistic. It is much weaker 
than the NN interaction, and has a small relativistic 
correction. 

In neutron matter at low densities the kaon energy 
deviates from its rest mass by the "Lenz" potential pro- 
portional to the scattering length 
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which is the optical potential obtained in the impulse ap- 
proximation. At high densities the kaon energy deviates 
from its rest mass by the Hartree potential: 



UJ H artre 



m K = p V K - n (r)d 3 r 
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-R 3 pV a , (7) 



which, as shown in H, is considerably less attractive. 

To estimate the transition from the low density Lenz 
potential to the high density Hartree potential we solve 
the Klein-Gordon equation for kaons in neutron matter 
in the Wigner-Seitz (WS) cell approximation. It allows 
an analytical calculations of the kaon energy which will 
be compared to numerical calculations on a cubic lat- 
tice. The cell boundary condition contains the important 
scale of the range of kaon-neutron correlations in matter, 
and gives the correct low density (Lenz) and high density 
(Hartree) limits. 

The RCKG equation for s-waves is 



/ tun + uj — V 
V rn N 



which has the solution 
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for r <R 
for r>R 



where n = ((cj + Vq) — m K )mN / {mN 



u) + Vo) and 

k 2 = {m 2 K — ui 2 )iriN I {rriN + uj). By matching the wave 
function u and its derivative at r — R we obtain 
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In the WS approximation 4>'(ro) = 0, where the cell 
size tq is given by the density p = (47rrQ/3) _1 . From Eq. 
(O) this implies 
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Eliminating the coefficient B/A from Eqs. ( ]To| ) and jTl| ) 

gives 



— tan(«;_R) 

K 



e 2k(R-r ) _ kr )/(l + fcr ) 
e 2fc(fl-r ) + (i _ kr )/(l + kr ) 



(12) 



which determines k and thus the kaon energy. The re- 
sulting kaon energy is shown in Fig. (Q). 

At low densities, r 3> R, the kaon energy is uj ~ 
and k ~ 0. Since also kr<j 1, we can expand the r.h.s. 
of Eq. and find 



tan(K_R) L, o l,. , 

^ '- = R+ -k 2 rl + -k 4 r 5 a + 

k 3 5 

Now we can extract the kaon energy 
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The linear part of Eq. (|TJ) is the Lenz potential, Eq. (|^), 
since uj + mx ~ 2m k at small p. The next order scales 
with p 4 / 3 and it becomes a quarter of the leading Lenz 
potential at p ~ po/16. This demonstrates that the Lenz 
potential of kaons is valid only at very small densities, 
and is of limited interest for kaon condensation. 

At the density where r — R equations (8) and (9) are 
solved by a constant <p implying uj = mx — Vo, the Hartree 
energy. At even higher densities the two-body potentials 
overlap with each other, and the Hartree approximation 
presumably becomes valid. It gives for r < R: 



m K - V Q 



(15) 



Note that this equation is valid in both relativistic and 
non-relativistic mean field limits; only the value of Vq is 
influenced by relativistic effects in the scattering process. 
The energies obtained with the Lenz and Hartree approx- 
imations are also shown in Fig. (Q). The cross-over from 
the Lenz to the Hartree limit takes place at rather small 
densities. 

If the kaon-nucleon potential has a short range repul- 
sive core of radius i? c , a stronger attractive potential — Vq 
is needed at R > r > R c to obtain the same scattering 
length. In Fig. ([!]) we show an example of the effect of 
a repulsive core. The repulsion is chosen such that the 
interaction has zero volume integral, i.e., the core po- 
tential is V c = V {R 3 /Rl - 1). We choose R = 1 fm 
and R c — R/2. In order to obtain the scattering length 
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ax-n — —0.41 fm, with non-relativistic kinematics, the 
attractive potential depth is Vq — 153MeV. At very low 
densities the kaon energy calculated with the WS approx- 
imation, follows the Lenz potential and is not affected by 
the presence of a repulsive core. At intermediate den- 
sities, tq ~ R, it is actually lower with than without a 
repulsive core. However, at higher densities the kaon en- 
ergy approaches uik , the Hartree limit for an interaction 
with zero volume integral. The presence of a repulsive 
core will thus further reduce the possibility of kaon con- 
densation. 
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FIG. 1. Kaon energy as function of neutron density. The 
results obtained with the WS method (Eq.12) are shown by 
full curve. At low densities they approach the Lenz result 
(Eq. (^), short-dashed curve) and at high densities they ap- 
proach the Hartree result (Eq.(jjj), long-dashed curves), while 
the dash-dotted curve gives the WS energies for the repul- 
sive core potential. The electron chemical potential p e (p) is 
shown (lower full curves) for the equation of state for beta 
stable nuclear matter ||[lo| without and with a transition to 
a mixed phase of quark matter for bag constants B — 122 and 
B = 200 MeV fm -3 . Quark matter admixtures reduce the p e 
as shown. 

In order to test the WS approximation we consider an 
extreme limit in which the nuclcons are infinitely massive 
and fixed in a simple cubic lattice. Using non-relativistic 
kinematics, it is simple to solve for the ground state of 
the kaon in this lattice. We evaluate the imaginary-time 
propagator exp[— Ht] for small r on a grid, and iterate 
^(t + At) = exp[— Hr^fr) until convergence to the 
ground state. This can be done efficiently by starting 
with a course grid and then using a finer mesh as the iter- 
ations proceed. The WS results for non-relativistic kine- 
matics and infinitely heavy nucleons are compared with 
the lattice results for a R = 0.7 fm square well potential 
with a,K-n = —0.41 fm. The structure of the lattice is 
not important; results for the bcc lattice, for example, 
fall between the WS and the exact results for the simple 



cubic lattice. In Fig. (Q) we show (u>k — m,k)po/ P which 
equals 2iraK- n Po/ m R m the Lenz and — Vo-R 3 47rpo/3 in 
the Hartree limits. We find little difference between the 
two sets of results throughout the range of densities con- 
sidered. 



-15.0 



> 

CD 



-20.0 



-25.0 



E 
I 

3 -30.0 



-35.0 



Hartree 








f / ^ - 

7/ 
// 


Wigner-Seitz 

• • Lattice 

— EE &«■„ = 1.54) 
- EE for, = 1.23) 




Lenz 


R = 0.7fm 



0.0 



2.0 
P'Po 



3.0 



4.0 



FIG. 2. Comparison of Wigner-Seitz energies, scaled by the 
factor po/p, with numerical exact results for a simple cubic 
crystal (dots), and Ericson-Ericson energies for £pro = 1.54 
and 1.23. These non-relativistic calculations use static, in- 
finitely massive nucleons. 

In Ref. Q] the kaon energy is corrected for correlations 
in the medium. This effect is analogous to the Ericson- 
Ericson (EE) correction for pions in a nuclear medium 
11 1 and the Lorentz-Lorenz effect in dielectric medium 
12 1 . It is interesting to compare the kaon energy loee 
obtained with the EE method for a cubic massive neutron 
lattice. It is given by M: 



luee ~ m K 



1 2nap 
vhk 1 — o£p 



3 C(r) 



(16) 



where C{r) is the nucleon-nucleon correlation function 
and ^/o is the inverse correlation length. In a simple cubic 
lattice £pTo = 1.54, while realistic neutron matter wave- 
functions calculated in || give £pr ~ 1.23 at densities 
> 2po- For comparison £ pro = 0.92 in neutron Fermi-gas. 

At lower densities there is little difference between the 
EE results for £pro = 1.54 and the WS or the exact re- 
sults. However, at higher densities, the luee is larger 
than the correct result. Note that the WS and the exact 
results depend upon the interaction range, while the loee 
depends only on the scattering length a. For larger values 
of i?, same a, the exact and WS energies are lower, fur- 
ther below the EE result, while for smaller R they would 
move up towards or even above the EE result. This is to 
be expected, since the EE approximation assumes that 
the kaon interacts with only one nucleon at a time, which 
is valid when ro 3> R- At high densities the EE method 
is not expected to be useful, however, it does seem to 
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be qualitatively good for R ~ 0.7 fm. The EE results 
are not excessively sensitive to the value of £p r o- In Fig. 
(0) we also show the results obtained with £pro — 1.23 
appropriate for realistic neutron matter. They cross the 
Hartree line at p ~ Apg. 

Self-energies in dilute systems can be expanded in 
terms of the scattering length (the so called Galitskii's 
integral equations 13|). Analogous results are obtained 
in chiral perturbation theory However, such expan- 
sions are valid only when the interparticle distance is 
much larger than the pair interaction range, so that only 
the scattering length matters. They generally do not 
have the correct high density limit, which depends upon 
the shape of the interaction in addition to the scattering 
length. 

The lowest order constrained variational method |l5| ] 
used to treat strong correlations in nuclear matter and 
liquid Helium is identical to the WS cell approximation 
employed here if the healing distance is chosen as r . 
These methods have the correct low and high density 
limits, and are meant to provide a good approximation 
over the entire density range. They are probably less 
accurate than the low-density expansions in the region 
where the expansions are valid. For example, consider 
the low density expansion of the WS energy (Eq. 14). 



To second order in the scattering length it corresponds to 
£pro = 1.8, much larger than the realistic values quoted 
above. Thus it is likely that at densities <C po the EE 
results for neutron matter shown in Fig. 2 for £,pr$ = 1.23 
are more accurate that the WS. However, the possible 
error in the WS results at p <C po seems to be < 5% by 
comparison. It is also possible to calculate corrections 
to the Hartree potential at high densities by coupling 
the kaon motion to phonons. Their estimate at p = 
4p is ~ —16 MeV for the square well potential with 
0.7 fm radius. This correction decreases as p increases. 
It is larger for Yukawa shaped potentials (same a and 
comparable radius) than the square well, however the 
Hartree potential is less attractive for the Yukawa shaped 
interaction. 

In conclusion, it seems that we can use the Hartree 
limit to estimate the kaon energy in matter at densities 
> 4po where kaon condensation may occur. On this ba- 
sis it appears from Fig. 1 that, for the electron potential 
p e {p) calculated from modern realistic iViV interactions 
without quark drops in matter, kaon condensation 
is unlikely up to p — 7po, which is just above the esti- 
mated range of densities possible in neutron stars. With 
the p e (p) used in Ref. |Q, which is larger, the condensa- 
tion could occur at p > 6po- If dense matter has ~ 10% 
protons, whose interaction with the K~ is believed to be 
twice as strong, the Hartree potential will be more at- 
tractive by ~ 10%. This together with the realistic p e , 
could reduce the condensation density to ~ 6.5po- If the 
K~N interaction has the Yukawa shape, or a repulsive 
core, that would push the condensation density higher, 
and if quark drops or hyperons reduce the value of p e (p) 
at higher densities as indicated in Fig.l, kaon condensa- 



tion becomes unlikely. 

We acknowledge Wolfram Weise for stimulating discus- 
sions. The work of VRP is partly supported by the US 
National Science Foundation via grant PHY 98-00978, 
the work of JC is supported by the U. S. Department of 
Energy under contract W-7405-ENG-36. 



[1 
[2. 
[3 
[4 

[5 
[6 
\-\ 
[» 
[9 
[10 

[11 

[12 
[13 

[14 

[is; 



D. B. Kaplan and A. E. Nelson, Phys. Lett. B291 (1986) 
57. 

G. Brown, C. Lee, M. Rho and V. Thorsson, Nucl. Phys. 
A 572 (1994) 693. 

V.R. Pandharipande, C.J. Pethick and V. Thorsson, 
Phys. Rev. Lett. 75 (1995) 4567. 

T. Waas, M. Rho and W. Weise, Nucl.Phys. A617, 449- 
463 (1997); T. Waas, N. Kaiser and W. Weise, Phys. Lett. 
B379 (1996) 34. 

R. B. Wiringa, V. Fiks and A. Fabrocini, Phys. Rev. C38 
(1988) 1010. 

R. Biittgen, K. Holinde, A. Muller-Groeling, J. Speth 

and P. Wyborny, Nucl. Phys. A506 (1990) 586. 

C. Lee, H. Jung, D. Min, M. Rho, Phys. Lett. B 326 

(1994) 14. 

N. Kaiser, P.B. Siegel and W. Weise, Nucl. Phys. A 594 

(1995) 325. 

A. Akmal, V. R. Pandharipande, & D. G. Ravenhall, 
Phys. Rev. C58 (1998) 1804. 

H. Heiselberg an d M. Hjorth- Jensen, Phys. Rep. (1999); 



nucl-th/9902033 



M. Ericson and T. E. O. Ericson, Annals of Phys. 36 
(1966) 323. 

G. Baym and G. E. Brown, Nucl. Phys. A 247 1975 395. 
S.T. Beliaev, Sov. Phys.-JETP, 7 (1958) 299. See also 
Fetter and Walecka, "Quantum Theory of Many Particle 
Physics", McCraw-Hill 1971. 

C. Lee, G. Brown, D-P. Min and M. Rho, Nucl. Phys. A 
585 (1995) 401. 

V.R. Pandharipande and H.A. Bethe, Phys. Rev. C7 
(1973) 1312. 



4 



